Quantum State Transfer in Spin-1 Chains 
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We study the transfer of quantum information through a Heisenberg spin-1 chain prepared in its 
ground state. We measure the efficiency of such a quantum channel via the fidelity of retrieving 
an arbitrarily prepared state and via the transfer of quantum entanglement. The Heisenberg spin-1 
chain has a very rich quantum phase diagram. We show that the phase boundaries are reflected 
in sharp variations of the transfer efficiency. In the vicinity of the border between the dimer and 
the ferromagnetic phase (in the conjectured spin-nematic region), we find strong indications for 
a qualitative change of the excitation spectrum. Moreover, we identify two regions of the phase 
diagram which give rise to particularly high transfer efficiency; the channel might be non-classical 
even for chains of arbitrary length, in contrast to spin- 1/2 chains. 

PACS numbers: 03.67.Hk,75.10. Jm,03.75.Mn 



Since the early works of Osborne et al. El and Oster- 
loh et al. , displaying an intertwined relation between 
entanglement and quantum phase transitions, the fields 
of condensed matter and quantum information have de- 
veloped a strong synergy, further motivated by the spec- 
tacular advances reached in the area of ultracold atomic 
physics and ion traps Among the broad scope of 
problems that nowadays can be addressed with ultra- 
cold atomic gases, spin models are particularly appealing. 
These relatively simple models exhibit the most funda- 
mental physics associated with magnetic ordering, criti- 
cality, and quantum phase transitions, but mostly lack an 
analytical solution. Though spin models have been con- 
structed as idealizations or toy models of real systems, 
ultracold atoms allow for an almost perfect realization 
of many of them. An example is the one-dimensional 
(ID) spin-1 system, which can be realized through con- 
fining an S = 1 spinor condensate in an optical lattice 
[3, HI- Restricting to nearest-neighbor interactions, the 
most general isotropic Hamiltonian for the spin-1 chain 
is the bilinear-biquadratic Hamiltonian (BBH) 
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Here Si = (Sf , Sf , Sf) are the spin operators on the ith 
site, and cos# (sin (9) gives the strength of the bilinear 
(biquadratic) coupling. The properties of the ground 
state as well as of the excitations are determined by the 
angle 8. The phase diagram is shown in Fig. [TJa). In 
the whole range — 37r/4 < 8 < x/2, the ground state 
is antiferromagnetic, i.e., has vanishing magnetization: 
M = (J^i Si) = 0. Since the Haldane conjecture that ID 
isotropic antiferromagnets with integer spin must have 
a unique massive, i.e., gapped, ground state with expo- 
nentially decaying correlations, the BBH has been exten- 
sively studied. The Haldane conjecture was rigorously 
proven for the AKLT point (tan 6* = 1/3), for which Af- 
fleck et al. explicitly constructed the ground state and 
proved the existence of a gap [f|. At 9 = ir/4 (Uimin- 



Lai-Sutherland point @, the system enters into a crit- 
ical (gapless) phase (ir/4 < 8 < tt/2) with unique ground 
state and diverging correlation length. At 6 = —ir/4 the 
gap vanishes Q, but it reopens for 9 < — 7r/4, where the 
systems enters into a dimerized phase. The exactly solv- 
able point 8 = — 37r/4 marks the border to the ferromag- 
netic phase (characterized by M ^ 0). The existence of 
a small spin-nematic region between the dimerized phase 
and 8 = — 37r/4 is actively discussed since a conjecture of 
Chubukov [lj|. 

In the field of quantum information, spin chains have 
been intensively studied regarding their usefulness as 
quantum channels EH Ell EH- Attention has been de- 
voted nearly exclusively to spin-1/2 chains, where it has 
been shown that a general quantum state can be trans- 
fered with relatively high fidelity between the two end- 
points of a ferromagnetic chain with nearest-neighbor in- 
teractions [111, E2l • As models with more complex ground 
states, spin-1/2 chains in the vicinity of a quantum phase 
transition EH, spin-1/2 ladders EH and Peierls distorted 
chains EH have been employed as quantum channels. 

Here we investigate the usefulness of the spin-1 chain 
as a quantum channel. To this aim, we study the trans- 
mission of an arbitrary quantum state, either pure or 
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Figure 1: (a) Quantum phases of the spin-1 chain, (b) 
Scheme for state transfer with an underlying ferromagnetic 
state, (c) with the channel initialized to the ground state of 
H(0), and (d) transfer of quantum entanglement. 
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entangled with a further subsystem, through the chain 
prepared in its ground state for the whole phase dia- 
gram. Compared to previous works studying either the 
ground state magnetic order or the excitation spectrum 
of the BBH, our approach combines both aspects simul- 
taneously. We find that the transfer does strongly de- 
pend on the nature of the ground state. When crossing 
the phase boundaries, the quality of the transfer changes 
sharply. Inside the Haldane phase (around the AKLT 
point) transfer is very inefficient. On the other hand, 
high transfer fidelities can be achieved by preparing the 
spin-1 chain in either the dimerized (gapped) or the crit- 
ical (gapless) phase. In those cases, transfer fidelities are 
clearly larger than for a ferromagnetic initial state. Fi- 
nally, we find a strong reduction of the transfer efficiency 
for 9 > — 37r/4, where a spin-nematic phase has been 
conjectured. This finding supports recent studies show- 
ing a qualitative change in the low-energy excitations as 
compared to the dimerized phase. 

Our scheme for quantum communication general- 
izes the one usually employed for spin- 1/2 systems 
(Fig. O^b)): we consider a chain of N sites with the 
first spin in an arbitrary state = Yl m =o ±i £ml TO )i 
(S z \m) = m\m), J2 m l£«i| 2 = 1) an d decoupled from the 
rest of the chain. The other A — 1 sites are, for a given 
9, prepared in the ground state |gSg)a...JV 01 H(9). The 
initial state reads \ip^(9,t = 0)) = <£> |gs e )2...jv, scc 
Fig. fflje). At t = 0, we abruptly switch on the interac- 
tion between the first and second spin and let the system 
evolve, obtaining \^{6,t)) = exp[-itH(9)]\ip i (9,0)) (h = 
1). At time t, the quality of the transfer of |£) to the last 
spin of the chain is evaluated by the fidelity (£\pn(&, 
of retrieving |£) at the end of the chain. The reduced 
state of site N is pjvj&t) = tri.. J v_i|V'f 
The channel fidelity [17| is obtained from averaging over 
pure states |£): 



F(9,t)= / dt(£\p N (0,t)\0, 



(2) 



where <i£ is the SU(3) invariant measure. For spin-1 al- 
ways 1/3 < F < 1. We define F{9) = F(9,t*), where 
t* is the time for which the perturbation arrives at the 
end of the chain for the first time {i.e., we ignore later 
maxima from multiple reflections on the boundaries). 

We also analyze entanglement transfer, considering 
that the spin at the first lattice site is entangled with 
a spin outside of the chain, say at lattice site (see 
Fig. [ljd)). As a particular case we take a singlet state 
|s)oi = (|l,-l)oi-|0,0) i + |-l,l)oi)/\/3, such that ini- 
tially \tp s (9,0)} = |s)oi<S>|gs )2...jv. As before, at t = Othe 
coupling between sites 1 and 2 is switched on, while site 
always remains uncoupled. We quantify the transfer of 
entanglement to the end of the chain by the logarithmic 
negativity [l|| for sites and N 



LN(M)=log 2 fc( 



(3) 



where /3 jv(M) = tri.,jy_i | tp s (9, t))( ip s (9, t) | , V denotes 
partial transposition, and \\A\\i = V trA^A. As defined 
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Figure 2: Channel fidelity F (open circles) and entanglement 
measured via logarithmic negativity LN (filled circles) for the 
transfer between the endpoints of a chain with 25 sites (ob- 
tained with MPS methods, using D = 20 ... 25). The dashed 
horizontal line denotes the maximal fidelity F c iass = 1/2 which 
can be obtained through classical communication. Vertical 
dashed lines indicate the quantum phase boundaries. (See 
Fig. [3]for a zoom into the region — 0.767T < 6 < —0.7n). 



here, < LN < log 2 (3) gives an upper bound for the 
number of spin-1/2 singlets that can be distilled from 
p 0N (9,t). Again we will use LN(0) = LN(0,i*). 

We have computed F(9) and LN(#)_for chains of up to 
N — 73 sites using MPS simulations Our results for 
25 sites and — tt < 9 < tt are shown in Fig. [2] The dif- 
ferent quantum phases of the model are well reflected in 
the transfer efficiency. Before discussing several interest- 
ing points in detail, it is useful to rewrite Eq. ([T]) in terms 

of two-site projectors P-j = Y^ m \&t , m) ij (St , m\ onto 
states with total spin St (to = St, — St) 



H{9) =JJ2 H l3 {9) =JJ2 *oP§ ] + XiPV + A 2 P, 



j(2) 



Here Ao = — 2cos0 + 4sin0, Ai = — cos6> + sin^, and 
A 2 = cos 9 + sin 9. We start our discussion at the border 
of the critical phase and follow the phase diagram in a 
counter-clockwise order. 

(a) Uimin- Lai- Sutherland point (9 = tt /A) and critical 
phase (V/4 < 9 < tt/2). Fidelity F and logarithmic 
negativity LN attain their maximum at 9 = tt/4, and 
decrease towards the ferromagnetic phase. For 9 = tt/4, 
the BBH takes a particular simple form 
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where W l} = (-1) 2S Est=o(~ 1 ) 5t ^ ( / t) is the P era " 
tor swapping sites i and j and ly is the identity op- 
erator. For N = 4, the initial state is | , 0{(vt/4, 0)) = 
(8 |gs T / 4 ) 2 34, with the trimer state 
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|t> 234 = -^=^(-1)^^^11,-1,0), (5) 
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Figure 3: Fidelity F (left column) and logarithmic negativity 
LN (right column) versus chain length N for (a,b) the spin-1 
chain at — tt/4; (c,d) the spin-1 chain at 8 — —tt/2; and 
(e,f) the spin- 1/2 Heisenberg chain. Dark circles correspond 
to the channel prepared in the ground state, open circles in 
(a,b) to a ferromagnetic initial state. F is fitted through an 
exponentially decaying function. Dashed lines indicate the 
classical communication limit -F c iass = 1/2 (-Fdass = 2/3) for 
spin-1 (1/2). (Results obtained with MPS simulations using 
up to D = 25 (D = 24) for 9 = tt/4 (-tt/2)). 



which has total spin zero [201 ] (P runs over permuta- 
tions of {2,3,4}, Pp permutes the sites accordingly). 
As \ip£(ir / 4, Jt = 71")) = |t) 123 ® |^)4, perfect transfer 
(F = 1, LN = log 2 (3)) occurs. For N > 4, trans- 
fer efficiency decreases (dark circles in Fig. [3] (a,b)). 
For small N the transfer is better for chains of length 
N = 1 mod 3, pointing to a trimerized order in the 
ground state. As N increases this difference vanishes. 
The average velocity v = N/t* at which the excita- 
tion propagates, grows as N increases. From simulations 
for N < 28, we extrapolate limjv^oo v w 1.59 Jsites, a 
value close to the velocity of sound in the infinite system 
v s = 2tt/ (3 V2) Jsites fx 1.48 Jsites [|. 

To emphasize the role of the initial state of the chan- 
nel, we use a ferromagnetic state \ipF(t = 0)) = (8> 
|1, 1, . . . , 1)2. ..n to compare with. As Hij swaps adja- 
cent sites, this reproduces the usual situation for spin-1/2 
chains [III, [13]. At time t 



N 



|^(7r/4,*))=€ 1 |l...l) + 5^7 ln (t) & 



(6) 



where |n,) represents the state with all spins in state 0, 
but the nth spin being in i = 0,-1. The corresponding 
probability amplitudes 7i n (i) can be calculated as an in- 
finite sum of Bessel functions [2l[. Already for N > 2, 
l7uv(t*)| < 1) an d thus also F < 1. As visible from the 
open circles in Fig. [3] (a, b), the transfer efficiency for the 
ferromagnetic initial state is much below the efficiency of 



the chain initialized to its ground state. 

(b) Ferromagnetic phase (9 G (7r/2,7r] U [—w, —3ir/4)). 
This region is characterized by A2 < A , Ai. The ground 
state has ferromagnetic order and broken rotational sym- 
metry: |gs fi (0,<p)) = ®£L 2 \U,<p)i- The state \l# llp )i 
has maximal spin projection in the direction specified by 

(p). Throughout the whole phase transfer efficiency 
is very small. The point 9 = tt/2 at the border to the 
critical phase (where for the finite chain the ground state 
is ferromagnetic) allows to identify two reasons: (i) for 
fixed (■#, tp), |£) = |0^ )V ) having vanishing z-projection in 
the corresponding direction is not transported; (ii) the 
transfer of |£) = | — is not via swaps (as for spin- 

1/2 or at 9 = Tr/A), but through an intermediate state. 

(c) Dimer phase (— Sir/A < 9 < —it/A). Inside this re- 
gion, F and LN increase strongly and reach a maximum 
at 6 = —n/2. At this point Ao < Ai = A2, thus the two- 
site ground state is a singlet. As (g> |s)23 ± |s)i2 <8> |C)3 
are both eigenstates of H(—ir/2) (with different eigenval- 
ues), for N = 3 perfect transfer occurs. For larger (odd) 
N the ground state of the last N — 1 spins is dimerized 
@, ^ e -i the expectation values of s^+i = |s),.i+i(s| 
are different on even and odd bonds. This is not the 
case for even A, and correspondingly transfer fidelities 
vary strongly between even and odd A, see Fig. [3jc,d). 
Note that H(—ir/2) = — H(%/2), i.e., the ground state of 
one model is the state with highest energy of the other, 
but fidelities observed in both cases are very different. 
This observation, which indeed can be made also at other 
points in the phase diagram, is a clear manifestation of 
the dependence of the transfer on the underlying mag- 
netic order and the nature of low excitations. 

Conjectured spin-nematic phase (—0.75ir < 9 < 
~0.72ir). Near the border to the ferromagnetic phase, 
Fig. [5] shows a remarkable drop-off in transfer efficiency. 
Fig. |4ja) shows a zoom into this region where a spin- 
nematic phase has been suggested [Til. This claim 
has been actively discussed recently but a complete 
characterization of this region is still under discussion 
@, |H [H, [H, [H, [H, [13. The dip for F and LN how- 
ever is consistent with the results of Lauchli et al. (26| 
and Porras et al. [27j , who found a qualitative change in 
the excitation spectrum for values — 37r/4 < 9 < 9c, with 
9c /tt w —0.7 . . .—0.67. To better characterize the ground 
state in this region, we calculate dimerization, nematic 
order parameter and magnetization (see Fig. H{b~d)). 

(d) Haldane phase (—tt/4: < 9 < tt/4). A large re- 
gion in this phase shows poor transfer with fidelities 
close to the lower limit F = 1/3. At the AKLT point 
(tan# = 1/3) this can be understood from the equiv- 
alence of the s pin- 1 chain to a model of two coupled 
spin-1/2 chains [22|: while the couplings between the two 
chains are completely symmetric, the AKLT ground state 
presents an asymmetry due to the boundary spin-1 /2s . 
This may lead to low transfer efficiency. 

Let us finally discuss quantitative features by compar- 
ing to the spin-1/2 chain without local engineering of 
couplings [2J|. Contrary to Ref. [ll|, where such a sys- 
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Figure 4: (a) F (open circles) and LN (filled circles) around 
the region for which a spin-nematic phase has been conjec- 
tured (for TV — 25). The transfer efficiency is strongly reduced 
for — 0.75-7T < 8 < 0.727T, as most clearly visible from LN = 0. 
To characterize the ground state we plot (b) the dimeriza- 
tion £>(jv-i)/2 (A = (Hi.i+i) - (Hi+i,i+2)), (c) the nematic 

order parameter Q = maxn Y^i=2 [((^n^i) 2 ) ~~ 2/3]/(AT — 1) 
(nn is the unit vector pointing in direction fi), and (d) the 
magnetization M = maxn "^f^i^Si) / (N — 1). 

tern was introduced, we always calculate F and LN for 
the first maximum, disregarding better values at (pos- 
sibly much) later times due to multiple reflections and 
constructive interference (3(j. F and LN for the spin- 
1/2 chain initialized to a ferromagnetic state and using 
a Heisenberg Hamiltonian are plotted in Fig. [3] (e,f). A 
good criterion to evaluate the efficiency of the channel is 
to compare its fidelity with the highest fidelity for trans- 
mission of a quantum state through a classical channel 
-Fciass = 2/(d+ 1), being d the dimension of the quantum 
system |17| . In particular, we can compute the maximal 
length of the channel, i.e., the maximal number of sites 



for which F > .Fciass- According to our numerical re- 
sults, for the spin-1/2 chain we obtain F < F c \ ass = 2/3 
for TV < 33. For spin-1 at 9 = —tt/2 we find that 
F > F c \ ass = 1/2 for all our simulations, i.e., for TV < 73. 
From fitting the fidelity through an exponentially decay- 
ing function, we extrapolate i 7 !^ = limjv— >oo F « 0.56, 
indicating that the channel is always superior to any clas- 
sical channel (we cannot exclude a different decay for 
larger chains, though TV = 73 is well above the dimcr 
coherence length TVq m 20 [H}). At = ir/4, fidelities 
are even larger. Considering chains of TV < 28 sites, we 
find that fidelity decays exponentially with asymptotic 
limit Fqo f=a 0.72 well above the classical value (again we 
cannot exclude a different decay for larger chains). 

Summarizing, we have studied state and entanglement 
transfer in spin-1 Heisenberg chains. We have shown 
that the quality of transfer, characterized by the aver- 
age fidelity and the logarithmic negativity, undergoes 
sharp changes at the phase boundaries. The critical 
(gapless) and the dimerized (gapped) phase have high 
fidelities; from extrapolating our data the channel might 
even be non-classical for any number of sites. In con- 
trast, transfer efficiency is significantly lower in the fer- 
romagnetic phase, and attains a minimum in the Haldanc 
phase (around the AKLT point) and in a small region at 
the border between the dimerized and the ferromagnetic 
phase, where the existence of a spin-nematic order is con- 
troversially discussed. 
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